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Abstract 

In this paper the first exact interior solution to Einstein's field 
equations for a static and non-uniform braneworld star with local and 
non-local bulk terms is presented. It is shown that the bulk Weyl 
scalar U(r) is always negative inside the stellar distribution, in con- 
sequence it reduces both the effective density and the effective pres- 
sure. It is found that the anisotropy generated by bulk gravity effect 
has an acceptable physical behaviour inside the distribution. Using a 
Reissner-Nordstrom-like exterior solution, the effects of bulk gravity 
on pressure and density are found through matching conditions. 
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1 Introduction 



The consequences of the braneworld theory [I] in general relativity have been 
studied with great interest during the past years j2] (see also a review paper 
[3] and references therein). Its possible effects on our observable 4D universe 
have been extensively studied through cosmological scenarios [I] (see also [5] 
and references therein) . The study of the consequences in astrophysics [5] has 
been mostly limited to the exterior region, even though it is well known that 
gravitational collapse could produce very high energies in the interior, where 
the braneworld corrections to general relativity would become significant [7] 
(see [8] for a recent study of black holes on the brane). 

In the astrophysics context, the general fact that a single 5D solution can 
generate different scenarios in 4D [9], makes the analysis in the exterior re- 
gion particularly attractive in searching for solutions beyond Schwarzschild's. 
Hence several scenarios which might be useful in predicting observable ef- 
fects from the extra dimension can be considered. On the other hand, the 
study on internal stellar structure in the braneworld remains unknown so 
far. The reason is very simple: the internal stellar structure is more complex 
that any other scenario, thus the high energy and nonlocal corrections to 
Einstein's field equation, which in general lead to a complicated indefinite 
system of equations in the brane, produce an even more complicated system 
of equations in the interior. Hence the construction and eventual study of 
internal stellar solutions in the brane is a difficult process, except when uni- 
form distributions are considered^. Consequently there are some important 
questions which remain without answer. For instance, the role played by 
density gradients as a source of Weyl stresses, and its eventual consequence 
on the gravitational collapse remain unknown. An internal consistent solu- 
tion would be useful in order to consider these issues. However, in general, 
finding a consistent solution in the brane, no matter which scenario is being 
considered, represents a challenge which final answer require more informa- 
tion on the 5D geometry and the way on how our 4D spacetime is embedded 
in the bulk. 

1 In the pioneer work of Germani and Maartens [lOj exact solutions for a uniform 
distribution were found, where it can be seen how much more difficult it would be to find 
a solution except for uniform stellar distributions. 
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There is a large number of studies on the consistency of braneworld mo- 
dels in the brane. For instance, the approach developed in the series of 
papers by Mukohyama [TT] has been extensively used in calculating gravity 
perturbations in a braneworld scenario. Early works on consistent linearized 
gravity [12], tests of consistency beyond linear order [13], stability and energy 
conservation [13], thus as low energy effective theory in the context of the 
braneworld [15]- [T7], have been extensively used as well. On the other hand, 
and from the point of view of a brane observer, there is an issue which must 
be faced in the search of any consistent solution in the brane; namely, the 
non locality and non closure of the braneworld equations [IE]- [22]. This is 
an open problem for which a solution requires a better understanding of 
the bulk geometry and proper boundary conditions^]. It is well known that 
the source of the non locality and non closure of the braneworld equations is 
directly related with the projection of the bulk Weyl tensor on the brane. 
Several ways have been taken to overcome this problem, most of them based 
in restrictions on the tensor S^. For instance a restriction which has proven 
to be useful consists in discarding the anisotropic stress associated to £^ u 
|24j . A different and more radical restriction on 8^ consist in imposing the 
constraint = 0. However this condition, which was initially used in some 
papers, is incompatible with the Bianchi identity in the brane [25]. 

In this paper there will be no direct restrictions on S^ u . Instead of this 
approach, the method developed in Ref. [26] shall be used to solve the non 
closure problem of braneworld equations. This method, based in the fact that 
any stellar solution on the brane must have the general relativity solution 
as a limit, is used in this paper to generate the first exact and physically 
acceptable interior solution for a nonuniform stellar distribution having local 
and non-local bulk terms on the brane. 

This paper is organized as follows. In Section 2 the Einstein field equa- 
tions and matching conditions in the brane for a spherically symmetric dis- 
tribution is reminded. In Section 3 a regular and physically acceptable exact 
internal solution is found. In Section 4 an analysis of the solution is carried 
out. In the last section the conclusions are presented. 

2 For a discussion about the role played by the boundary conditions and the gravitational 
instability on the brane see [23] 
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2 The field equations and matching condi- 
tions 



The Einstein field equations on the brane may be written as a modification 
of the standard field equations [T8]'[T9] 
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where A is the cosmological constant on the brane. The energy-momentum 
tensor has new terms carrying bulk effects onto the brane: 
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here a is the brane tension. The new terms S^ v and £^ u are the high-energy 
corrections and KK corrections respectively, and are given by 
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being U the bulk Weyl scalar and and the anisotropic stress and 
energy flux respectively. 



We consider a spherically symmetric static distribution, hence = 



and 
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where r M is a unit radial vector and = — u^u u the projection tensor 
with 4- velocity u^. The line element is given in Schwarzschild-like coordinates 
by 

ds 2 = e u dt 2 - e x dr 2 - r 2 (d6 2 + sin 2 6d<f) 2 ) (6) 
where v and A are functions of r. 

The metric has to satisfy ([!]). In our case with A = we have: 
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Pi = -y(p + p), (10) 

where /i = df/dr and fc 2 = 87r. The general relativity is regained when 
a^ 1 — > and ( TTOi) becomes a linear combination of (J7j)-([9l). 

The Israel-Darmois matching conditions at the stellar surface S give 

[G>1s = (11) 

where = /(r) |#+ —f(r) \r- Using (TTTi) and the field equation (TjQ) with 
A = we have 
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where = /(r) \ r= R. The equation (fH]) gives the general matching condi- 
tion for any static spherical braneworld staio HOI ■ When a^ 1 — ► we obtain 
the well known matching condition pr = 0. In the particular case of the 
Schwarzschild exterior solution U + = V + = 0, the matching condition (TTJ 
becomes: 
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3 The general matching conditions on the brane for a spherically symmetric vacuum 
region embedded into a cosmological environment can be seen in |20j 
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Thus the matching conditions do not have a unique solution on the brane. 
It is easily seen that the field equations (J7])-([HD can be written as 
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The equation ([16]) actually represents an integral differential equation for 
the geometrical function A(r), something completely different from the gen- 
eral relativistic case, and a direct consequence of the non locality of the 
braneworld equations. The only solution known for this equation is given as 
[26] 
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The function H(p,p,v) measures the anisotropic effects due to bulk conse- 
quences on p, p and v, and its physical meaning will be useful in the searching 
of exact solutions. This will be addressed in the next section. 
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3 An exact solution 



So far we have the interior Weyl functions V and U plus the geometric func- 
tion A(r), respectively given by (iTTj) . (|T8|) and (I2TI) . and three unknown 
functions {p(r), p(r), ^(r)} satisfying one equation, namely, the conservation 
equation (flOl) . Therefore it is necessary to prescribe additional information 
to close the system. First of all, to ensure the correct limit at low energies, 
the following constraint is imposed on the brane 



The constraint (1251) has been proven to be useful in finding solutions which 
possesses general relativity as a limit [26J, and has a clear physical interpre- 
tation: eventual bulk consequences on p, p and v do not produce anisotropic 
effects on the brane. This constraint ensures that a braneworld solution be 
consistent with general relativity, as is shown in Ref. [26]. Unfortunately this 
constraint is not enough, and one additional condition must be imposed. 
Therefore, from the point of view of a brane observer, many solutions are 
possible. However not all of these solutions are of physical interest. Hence 
the brane observer has to impose a condition in the brane which must lead to 
a physically acceptable solution, namely, regular at the origin, pressure and 
density defined positive, well defined mass and radius, monotonic decrease 
of the density and pressure with increasing radius, etc. All these conditions 
reduce enormously the possible solutions (p, p, u) to ( TTOT) and ( 1251) . even more 
in the searching of an exact solution. 

It might be logical to think of the Schwarzschild condition e u = e _A to 
close the system. However this would produce a very complicated integral 
differential equation for A, as can be seen through (12T1) . On the other hand, 
it is clear that to find an exact expression for the geometric function A(r), 
given by fl2Tj) . an analytic solution for fl23|) is needed. Keeping this in mind, 
a huge and simple family of exact solutions for ( |23i) is considered, given by 



which is characterized by the constants A, C, m and n. The constants A 
and C are eventually found through matching conditions, whereas m and n 
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are parameters to be used in the searching of an exact solution for p(r) and 
p(r). 

Using (|2"oT) in (FlUi) the pressure is found in terms of the density 



2B - VACmn J ' r m ~ l (1 + C V™)? -1 p ( r ) dr 
~ 2y/A(l + Cr m )^ ' 



with B a constant of integration. 

Now using (126]) and (|27l) in the constraint (|25l) . the following integral 
equation for the density is obtained 
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m\ 2 



V^4(l + Cr 

4 7T (2 + C (2 + m n) r m ) (/ r r 2 p(r) rfr - r 3 p(r)) _ 2 = 
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Hence p(r) can be found by (!28|) and then p(r) through (|2Tj) . In this stage 
it is important to stress that the goal of this paper is to find a simple exact 
braneworld solution, hence we are not interested in the most general solution 
p(r) to the integral equation (1281) . Indeed, a particularly useful solution for 
(|27|) and ( |28l) . leading to an exact braneworld solution through ( |2TI) . will be 
constructed. This is shown below. 

The first step is to examine the complicated bulk contribution to A(r) in 
(pTj) . When the constraint (|25|) is imposed to ensure general relativity as a 
limit, the geometric function A(r) underwent a great simplification, leaving 
only high energy corrective terms. Then when ( |26l) is used, the complicated 
integral expression in fl2TT) is reduced even more, leading to 

IT pj / „ s r r r. i 4+n(m-3) 

/ 2- ( p 2 + 3pp ) dr = / 2r 2 (l + Cr m ) n " 1 [4 + C(4 + mn)r m ] *+™- 

(p 2 + 3pp) dr. (29) 
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Finding a physically acceptable exact solution for (1271) and (J2HD, which can 
produce an exact solution for (1291) . represents a difficult process. Hence the 
integral expression (1291) has to be simplified as much as possible. Indeed, the 
integral ( 1291) is tremendously simplified when 4 + n (m — 3) = 0, leading to 

T 7irr~~ — 2T (p 2 + 3pp) dr = f 2r 2 (1 + C r m )^ (p 2 + 3pp) dr. (30) 

The next step will be to consider in ( 1271) a simple ansatz for p(r) which is 
capable of producing a simple expression for the pressure, and then to use 
( 1281) to fix the parameters of the ansatz. The idea is to obtain a solution 
for both p(r) and p(r) as simple as possible such that the integral (1301) has 
an analytic expression. If this is accomplished, the geometric function A(r), 
given through (121]) . and both interior Weyl functions V and U, given by ( 1T71) 
and ( 1T81) respectively, will have exact expressions. 

By simple inspection of ( 1271) it is easy to see that a convenient ansatz for 
p{r) can be written by 

p(r) = (1 + Cr m )"W 2+1 ' ]T a s r s . (31) 

On the other hand, since (1261) must be regular at the origin r = 0, m has 
to be positive. The constant n, which is given by n = is considered 

positive to obtain a pressure with a physically acceptable behaviour. Thus 
in our case m satisfies < m < 3. Taking m = 2 and using (13T|) in (l27j) . it 
is found that p(r) has a relatively simple expression free of special functions 
when @ 

(a + a 2 r 2 + a 4 r 4 ) 
= (i + Cr-f ■ (32) 

Hence the following expression for the pressure is found 

(r) = C + Cir 2 + C 2 r 4 + 2 (2a 4 -a 2 C) (1+Cr 2 ) log(l + Cr 2 ) 



4 Using an algebraic manipulator, it is not complicated to realize that keeping a generic 
value for m and s produce a solution to (|27| having special functions. This special functions 



eventually would make impossible an exact solution to A(r) through ([2T 
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where 



C 2 = -2a 4 C 2 . (34) 

The expression ( 1331) will hardly produce an exact function to (!30|) unless the 
logarithmic function be removed. Thus, 

«4 = 7^2- (35) 

Now using (132p . the constraint ( 1281) is written as 



v 7 ^ ^-24C 2 + 66a 2 7r + 36a C7r+ ^^j^ ) r + v 7 ^ (-20 C 3 + 70 a 2 C i 
- (9 a 2 - 2 a C) 7T (5 + 9 C r 2 ) arctan(v / C r) = 0, 



hence it is found that 

9C 2C 2 C 3 UCVI 
a = -z—; a 2 = — — ; a 4 = — ; £ = . (37) 

7 7T 77T 77T 7/T 

Using ( 1341) and ( 1371) in (}3"2l) and ( |33|) . a simple and physically acceptable 
expression is found for both the pressure and density. Thus the solution for 
(TTU]) and ([25]) is finally written as 1 

e u = A(l + Cr 2 f (38) 



= C(9 + 2C l!± ^) 

V 7 7tt(1 + Cr 2 ) 3 V 7 

and 

, N 2C(2 - 7CV 2 - CV) . , 

p(r) = 7,(1 Tm - • (40) 

leaving A and C to be determined by matching conditions. 



5 The expression for v and p were found in [29] for a perfect fluid in the context of 
general relativity. However the solution shown there is not physically acceptable due to 

Gr _L 
r T 
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Using (135|) - (PiO|) in (12"TT) a regular and well defined solution for A(r) is 
obtained 

e-™ = 1 - (41) 
r 

where the interior mass function rh is given by 
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with m(r) being the general relativity interior mass function, given by the 
standard form 
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hence the total general relativity mass is obtained 
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where R is the radius of the distribution. 

Using (fTTjl and ({TBI a regular solution for the interior Weyl functions is 
obtained 

32 

441r 3 (l + CV 2 ) 6 (l + 3CV 2 ) 2 

+ 16533C 3 r 6 + 12660C 4 r 8 + U6C 5 r 10 - 120CV 2 + 9C 7 r 14 ' 
-60^(1 + CV 2 ) 3 (3 + 26CV 2 + Q3C 2 r 4 )arctg(VCr)} , (45) 
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U(r) 



32 
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+6186C 3 r 6 - 373C 4 r 8 - 219CV - 18C 6 r 12 
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(46) 



The expressions fl38l- f|40l with f|45|) and ()4"6l) represent an exact analytic 
solution to the system (J71)- (JT0|) . 
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4 Analysis of the solution 



As can be seen by figure 1, the scalar function U{r) is always negative inside 
the stellar distribution, with a maximum negative value at the origin r = 0. 
This situation may be explained through the general expression for U(r), 
given by Eq. (|T8|) . It shows two "sources" for U(r): the first kind given by 
the first two terms in the second hand side of Eq. ({TBI) , which are high energy 
corrections, always negative. The second kind is given by the remaining 
terms, which clearly represent an anisotropic expression, which is always 
positive. Hence when the anisotropy projected onto the brane is not high 
enough, the dominant high energy terms produce a negative scalar function 
U(r), which is the case presented here. This negative scalar function reduces 
both the effective density and the effective pressure, as can be seen by the field 
equations (jT])-®- On the other hand, the anisotropy inside the braneworld 
star is shown by the figure 2. It increases until reaches a maximum value, then 
decreases until V = at r = 0. This is directly connected with the correction 
for A proportional to high energy terms shown in ( l2Ti) . This correction is the 
only bulk effect underwent by the metric when the constraint (1251) is imposed, 
therefore it represents the only source for V, as can be clearly seen through 

<m- 

The bulk contribution to p, p and v is found by matching conditions, 
where the assumption of vanishing pressure at the surface will be dropped 
[27j'[28j. As Schwarzschild is not the only possible static exterior solu- 
tion, we have many scenarios to consider. For instance let us consider the 
Schwarzschild exterior solution 

e v+ = e~ x+ = 1 - — ; U + = V + = 0. (47) 
r 

The matching condition [ds 2 ]^ = at the stellar surface £ yields 

A = (l-^)(l+CR 2 )-\ (48) 
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(49) 



(1 + CR 2 ) 2 (1 + 3CR 2 )VCR\ ' 
and using (|15[) it is found that C must satisfy the condition 

irR [168CR 2 + 252{CR 2 ) 2 - 18A8{CR 2 ) 3 - 4032(Ci? 2 ) 4 - 2352(Ci? 2 ) 5 - 252(C7? 2 ) 6 ] 
+±C7? [240 + 2749C7? 2 + 5276(C7? 2 ) 2 - 266(C7? 2 ) 3 - 372(OR 2 ) 4 - 27(OR 2 ) 5 ] 

-±240\/Carctg(VCR) (1 + 11C7? 2 + 19C 2 i? 4 + 9C 3 i? 6 ) = 0, (50) 

hence solving (l5T?j) C is found as a function of the brane tension a. 

In order to find the bulk contribution to p and p we need to find C 
satisfying ( 1501 . To carry out this the following solution is proposed 

C = C + 5, (51) 

where Cq is the general relativity value of C, given by 

°° ~ ~2R^' (52) 

which is found using the condition p(R) = in (]40l . In this sense S represents 
the " bulk perturbation" of the general relativity value of C. Using (15ip in 
fl50|) we have at first order in a -1 
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(53) 



p(Co + S)=p(C )+5^\ c=Co , (54) 



The pressure can thus be found expanding p(C) around C ( 

. dp 

which leads to 



_ 2C (2 - 7C r 2 - C 2 r 4 ) 4 (1 - 9C r 2 + 2C V) 

Mrj "77 (1 + C r 2 ) 3 + 7^ (1 + C r 2 ) 4 d ' (55j 

By the same way the density is found to be 

. Co (9 + 2C r 2 + Cy) 1 ( 9 -14C r 2 + C 2 r 4 ) r 

P(r) = 77 (1 + C r 2 )3 +7^ (1 + C r 2 ) 4 * (56) 
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However at the surface and for any arbitrary R always we have 



4 (1 - 9C R 2 + 2C 2 R*) ^ 
P{R) = T, (1 + Q^ 2 ) 4 6 < °' 



(57) 



Hence the Schwarzschild exterior solution is incompatible with the interior 
solution found here. Thus a different exterior solution must be considered. 



Using now the Reissner-Nordstrom-like solution given in [30] 
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and considering the matching condition [ds 2 ]s = at the stellar surface E, 
we have 

A(l + CR 2 ) 4 = l-^ + ^, (60) 
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(1 + CR?f (1 + 3CR 2 ) VCR 
and using (fl4l) we obtain 

9 = 147(1 + C^f(l + 3CJP) ^ + ^ + 22C ^ 4 + 3C3fi6) 

84nR 2 (l + Ci? 2 ) 2 + -(-240 - 2749C7? 2 - 5276C 2 i? 4 + 266C 3 J R 6 
a 

+372C 4 i? 8 + 27C 5 i? 10 )) + -240VC(1 + C7? 2 ) 2 (l + 9CR 2 )arctg(VCR) 

' a 



(62) 



The constants M. and q are given in terms of C through equations (l6"Tj) and 
(|62|) respectively, and C may be determined by ([601 if A is kept as a free 
parameter, which can be used to find a physically acceptable model. However 
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we have to be aware of the fact that A has a well defined general relativity 
value, named Aq, which is given by d60j) at a' 1 = 

A (1 + C R 2 ) 4 = 1-™. (63) 

In this sense the free parameter associate to A, which will be used to obtain 
an acceptable model, is the "bulk perturbation" of A given through 

A = A + e. (64) 



Using ([511 and (I64|) in (|60|) we obtain 



{A + e)[l + (Co + 5)R 



2i-i 



2M 
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q_ 

R 2 ' 



(65) 



Evaluating the expressions (IBTT) and (|62l) at C = Co + 5 and keeping li- 
neal terms in a~ l , the equation (j65p leads to the explicit form to the bulk 
perturbation of C, which is written as 
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with 



a [a 



a 



2\ 2 Co 

7/ 7T 

80 



240 + 589C i? 2 



25C2.R 4 



41C^ 6 
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(67) 



Hence giving by hand the perturbation e underwent by A due to the extra 
dimension, it is possible to obtain 5 and thus the bulk consequences on p and 
p through fl55|) and fl56l) . The figure 3 shows the behaviour of the pressure 
in both the general relativity and braneworld case. It can be seen that the 
bulk gravity effect reduces the pressure deep inside the distribution, but the 
situation changes for the exterior layers, where the matching conditions lead 
to p 7^ at the surface. 
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U(r) 




Figure 1: The scalar function Uir) for a distribution with R = 5. Uir) is 
always negative in the interior, hence it reduces both the effective density 
and effective pressure. 
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Figure 2: Behaviour of the anisotropy Vir) inside the stellar distribution 
with R = 5. 
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Figure 3: Qualitative comparison of the pressure p(r) in general relativity ( 
p(R) = ) and the braneworld model ( p(R) ^ ) with R = 5. 

5 Conclusions 

In the context of the braneworld, a spherically symmetric, static and non- 
uniform stellar distribution with Weyl stresses was studied. The method 
developed in Ref. [26J, which is based in the fact that any stellar solution 
on the brane must have the general relativity solution as a limit, was used 
to overcome the non locality and non closure of the braneworld equations. 
Hence there was no direct restriction on the projected £^ v Weyl tensor on 
the brane, which is a method usually used. 

By prescribing the temporal metric component goo , the first exact and 
physically acceptable interior solution to Einstein's field equations for a static 
and non-uniform braneworld star was found. It was shown that this solu- 
tion is incompatible with the Schwarzschild's exterior metric. Using the 
Reissner-Nordstrom-like solution given in Ref. [30] , the effects of bulk grav- 
ity on pressure and density were found through matching conditions, where 
the assumption of vanishing pressure at the stellar surface was dropped. It 
was found that the bulk gravity effect reduces the pressure deep inside the 
distribution, but the situation changes for the exterior layers as a direct con- 
sequence of matching conditions, in agreement with the previous study in 
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Ref. 

It was found that the Weyl scalar function U(r) is always negative inside 
the stellar distribution. In consequence it reduces both the effective density 
and the effective pressure. On the other hand, the anisotropy inside the 
braneworld star, which is directly connected with the deformation underwent 
by A due to bulk gravity effects, has an acceptable physical behaviour. 

The exact solution found in this paper was possible as a direct conse- 
quence of the constraint H(p, p, v) = 0. This essentially allows us to simplify 
the solution for the geometric function A(r) by eliminating some anisotropic 
effects on the brane. As can be seen by Eq. (fTTl) . the source of V is the 
deformation underwent by the geometric functions A(r) and v(r) due to bulk 
consequences on the brane. However when the constraint H(p, p,u) = is 
imposed, the bulk effect on z/(r) does not produce any anisotropic conse- 
quence, leaving thus the deformation underwent by A(r) as the only source 
of anisotropy on the brane. Furthermore, this constraint reduces the defor- 
mation of A(r), leaving a corrective term proportional to high-energy effects 
of bulk gravity, as can be clearly seen through the equation (T2T1) . Since the 
constraint H(p, p,v) = removes all possible sources of anisotropy except 
for the corrective term of A(r) proportional to high-energy effects, it follows 
that the anisotropic effect of the bulk on the brane is reduced to its minimal 
expression, hence the constraint imposed represents a condition of minimal 
anisotropy on the brane. 

The condition of minimal anisotropy, represented by the constraint H (p, p, z/)=0, 
is not only a direct path to avoid the loss of the general relativity limit, but 
also a natural way to reduce the degrees of freedom on the brane. This is 
an enormous simplification which has proven to be useful in searching physi- 
cally relevant exact solution on braneworld. Furthermore, this method works 
equally for analytic as well for numerical methods. Hence the condition of 
minimal anisotropy might help in the search of not only analytic physically 
acceptable, but also numerical solutions when a non-uniform distribution is 
considered. Thus the role played by the density gradients as a source of Weyl 
stresses in the interior could be studied. 

The work developed in this paper represents the point of view of a brane 
observer. Hence the solution found here, a physically acceptable one, does 
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not ensure that the bulk eventually constructed will not be plagued with 
singularities or any other problem. However, since the condition of minimal 
anisotropy on the brane ensures the correct limit at low energies, it could be 
used when the bulk configuration is investigated, thus some general features 
of the five dimensional bulk might be elucidated 0. This is currently been 
investigated. 
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